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1. Introduction 

Supersymmetry (SUSY) requires the existence of two scalar partners (squarks) for every quark 
corresponding to its two helicity states. Many models of SUSY predict that, due to large mixing, 
the lightest squark could correspond to one of the mass eigenstates of the third generation with 
niq < 500 GeV, thus allowing for the production of stop pairs at a future Linear Collider 
operating below 1 TeV. A “threshold scan” of the total cross section line-shape at such facility 
will yield precise measurements of the stop mass, lifetime and couplings, in close analogy to the 
program carried out in threshold studies for the top-anti top threshold [|^. It is then clear that 
an accurate description of the shape of the squark production cross section at energies close to 
threshold will be mandatory to put on equal footing the theoretical and experimental uncertainties. 

Close to threshold squark pairs are produced with small velocities in the center of mass frame, 
mv^ = — 2m, m being the heavy squark mass. The low-energy QCD dynamics of squarks is 

based on standard QCD if we assume that the gluino is not much lighter that the electroweak scale. 
The multi-gluon exchange of time-like gluons produces singular terms at threshold oc (a^/v)” and 
oc (a^lnv)” at the n-loop order in QCD perturbation theory, which forces one to carry out a double 
expansion in and v and count v ~ to define the power counting in this regime. The sin¬ 
gular terms mentioned above must be thus termed as leading order contributions and need to be 
summed up to all orders in For top-antitop pair production at threshold, vNRQCD ("velocity 
NRQCD”) [§, ^ has been shown to be the proper framework to consistently perform the summa¬ 
tion to NNLL order [Q], which means to account up to a^v, v^) corrections to the leading- 

logarithmic (LL) order oc L^:(o;i/v)^^,(ailnv)'. An analogous effective field fheory describing fhe 
non-relafivisfic inferacfion befween pairs of colored scalars can provide fhe needed ingredienfs for 
a renormalization group improved compufafion of squark pair production af fhreshold. A complefe 
summation af NLL order requires fhe knowledge of fhe Wilson coefficienls of pofenfials af LL or¬ 
der and of nonrelafivisfic producfion currenfs af NLL order ii- The effecfive Lagrangian for fhe 
scalar version of vNRQCD (fhaf will be called simply vNRQCD in fhe following) and fhe relevanf 
mafching condifions and anomalous dimensions are ouflined in fhis leffer. We refer fhe reader fo a 
recenf work ||^ for a more defailed analysis and explicif formulae. 

2. Effective theory set-up 

The vNRQCD effecfive fheory is formulafed by including only fhose quark and gluonic de¬ 
grees of freedom which can become on-shell for energies below m. The on-shell degrees of freedom 
correspond fo gluons {Aq,A^) and ghosfs (cq, c) wifh sofl ~ (mv,mv) and ullrasofl ~ {mv^,mv^), 
energy and momenfa, and pofenfial heavy squarks (t/Tp) and anfisquarks (Xp), wifh energy ~ mv^ 
and soff fhree momenfum p ~ mv. Massless quarks and squarks can also be included wifh soff and 
ullrasofl componenls in vNRQCD. The dependences on sofl momenfa of fhe heavy squark and sofl 
gluons appear as labels on fhe fields and an explicif coordinate dependence refers only fo ullrasofl 
flucluafions 0. All off-shell effecls such as fhose from hard squarks and gluons, potential gluons, 
and sofl quarks are accounled for by on-shell mafching of vNRQCD fo full QCD af fhe hard scale 
m. The basic vNRQCD Lagrangian can be separaled info ullrasofl, sofl and pofenfial componenls 
+ + The ullrasofl piece of fhe effecfive Lagrangian describes fhe inleraclions of 
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Figure 1: Effective theory diagrams for the interaction of a heavy squark with an ultrasoft gluon. 


ultrasoft gluons with squarks (Fig. |I|)a, which are multipole expanded to separate ultrasoft mo¬ 
menta from larger momenta, and also contains the kinetic energy contributions for the squarks and 
the ultrasoft fields: 
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where G^'' is the ultrasoft field strength tensor and the covariant derivative, + ijifjguA^, 

only involves the ultrasoft gluon field and the ultrasoft gauge coupling = gu{}^u)- The factor 
of the ultrasoft renormalization scale Hu that appears in the covariant derivative is determined by 
the requirement that the kinetic terms in the vNRQCD action are of order v** ^ in r/ = 4 — 2e 

dimensions, and we employ the same power counting formalism as in the heavy quark case li|l- 
Likewise, the velocity scaling determines the factors of the soft renormalization scale jXs accom¬ 
panying the soft gauge coupling gs = gs{lJ-s) in the soft and potential interactions (Figs. ||b,c). The 
existence of two different subtraction scales guarantees that each of loop in the EFT receives the 
proper renormalization scale according to the three-momenta flowing through it. A crucial fea¬ 
ture of vNRQCD is that fis and Hu are required to be correlated according to the non-relativistic 
energy-momentum relation of a heavy particle, i.e. jXu = jj.^/m = mv^, where v is the dimension¬ 
less renormalization parameter of vNRQCD. Lowering v from the hard matching scale (v = 1) 
to a scale of order the squark velocity sums all logarithms of the soft and ultrasoft scales into the 
Wilson coefficients of the vNRQCD operators, and renders vNRQCD matrix elements small. 

The vNRQCD Lagrangian also contain operators that describe interactions of quarks with 
soft gluons (Fig- lilb), that arise when soft quarks are integrated out, and potential-like squark- 
antisquark interactions fifp (Fig. |^), originating from potential gluons and other off-shell modes. 
Setting aside spin-dependent contributions that are absent for scalars, the matching coefficients of 
fif’s in the scalar theory agree with the corresponding ones in the quark theory due to the equiva¬ 
lence of the heavy quark and heavy scalar effective Lagrangians (HQET and HSET) at order 1/m. 
Differences arise at NNLO because the Darwin term vanishes at the hard scale in HSET [||]. The 
coefficients in fif’s are only renormalized due to soft interactions [^, and the structure of the 
resulting soft divergences is identical to the one in HSET. The potential Lagrangian has terms 


= -mI""L V (p,p') YpVpX*p'X p + ■■■ 

n n' 


y(p,p') = (r^®f^) 
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The coefficients iCFT) obtained by matching to the full theory Born scattering diagrams (e.g. 
for the Coulomb potential shown above, yf Hi) = 47rai(m)), and they differ in general from the 
corresponding results in the quark case. The EL renormalization of the potentials are determined 
by the effective theory one-loop graphs shown in Eig. Since to the order we are interested in 
interactions with ultrasoft gluons are not sensitive to spin, the 1-loop ultrasoft renormalization is 
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Figure 2: (a) Diagrams for the computation of the LL anomalous dimension of potentials, (b) NLL running 
of the normalized Wilson coefficient for S-wave and P-wave currents for m = 220 GeV (solid and dashed 
lines) and for m — 500 GeV (dot-dashed and dotted lines). 


the same in vNRQCD for both quarks and squarks. There are a number of operators in the potential 
and soft parts of vNRQCD with vanishing matching conditions at V = 1 which however become 
non-zero for v < 1 due to ultrasoft renormalization [||]. 


3. Anomalous dimensions for production currents 

At leading order in the non-relativistic expansion heavy squark-antisquarks pairs are produced 
in S- and P-wave states in 77 and collisions by the currents /s p = t/r* and Jpp = V^p p;^*p 
respectively, with corresponding Wilson coefficients Ci(v) and cp{v). The total cross sections 
77 —> qq) can be written in terms of the imaginary part of the correlators built from the 
non-relativistic currents and the products of their Wilson coefficients. The NLL running of the Wil¬ 
son coefficients is determined by the two-loop anomalous dimensions of the currents, as calculated 
in Ref. [^. The results are displayed in Fig.^. The scale variation for the S-wave coefficient is 
found to be much stronger than for the P-wave coefficient and the maxima slightly decrease with 
the heavy scalar mass and also move towards smaller values of v. The remaining scale depen¬ 
dence in the NLL cross section comes from the non-relativistic correlators that are proportional to 
the zero-distance Green’s functions obtained by solving the Schrodinger equation with the NLL 
improved Coulomb potential. 
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